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Abstract. In this paper we discuss local averages of the energy density for the non-minimally 
coupled scalar quantum field, extending a previous investigation of the classical field. By an 
explicit example, we show that such averages are unbounded from below on the class of Hadamard 
states. This contrasts with the minimally coupled field, which obeys a state-independent lower 
bound known as a Quantum Energy Inequality (QEI). Nonetheless, we derive a generalised QEI 
for the non-minimally coupled scalar field, in which the lower bound is permitted to be state- 
dependent. This result applies to general globally hyperbolic curved spacetimes for coupling 
constants in the range < £ < 1/4. We analyse the state-dependence of our QEI in four- 
dimensional Minkowski space and show that it is a non-trivial restriction on the averaged energy 
density in the sense that the lower bound is of lower order, in energetic terms, than the averaged 
energy density itself. 

PACS Numbers 04.62. +v 



1 Introduction 

For more than 30 years it has been known that the stress-energy tensor of the classical 
scalar field, obtained from the Lagrangean 

L[<P\ = l -{V<P?- l -{m* + iR)<p\ (1) 

does not satisfy the weak energy condition (WEC) at non- minimal coupling, i.e., £ ^ 
(see [1] and [2] for a simple example in Minkowski space). Naturally, this raises the question 
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of whether there are any restrictions on the extent of WEC violation and whether this field 
could support exotic phenomena (e.g., violations of the second law of thermodynamics) 
that depend on macroscopic spacetime regions of negative energy density. 

We showed in [3] that under certain conditions one can find lower bounds for local 
averages of the energy density of the classical non-minimally coupled scalar field. This 
paper extends our analysis to the case of the quantised field. 

As is well known, quantum field theories obeying the Wightman axioms necessarily 
violate the WEC [I] and in this respect, the non-minimally coupled scalar field resembles 
the situation at minimal coupling. However, as we show, there are differences. For example, 
consider the case where £ > 0. Given any bounded subset O of Minkowski space and an 
arbitrary constant p > we w ih construct a Hadamard state \1/ in which the expected 
energy density (p)^ is less than — p throughout O. Therefore, non-trivial local averages 
of the energy density, (p(f))*, where f is a non-negative test function, are unbounded from 
below on the class of Hadamard states. 

By contrast, expectation values of the averaged energy density of the minimally coupled 
scalar field are bounded from below on the class of Hadamard states [5] . The latter bound, 
known as a quantum energy inequality (QEI), can be written in the form 

<p(f)>* > -0(f) (2) 

for all Hadamard states where 0(f) is a constant (see [HI [3, [8] for reviews and further 
references concerning such QEIs). Clearly, the non-minimally coupled field cannot satisfy 
a QEI of this type. However, we will show that it obeys a generalised QEI of the form 

<p(f)>* > -(0(f))* (3) 

for all Hadamard states where 0(f) is now allowed to be an unbounded operator, which 
turns out to involve the Wick square of the field. This bound will be proved for averaging 
along time-like geodesies in general globally hyperbolic spacetimes. Precise statements and 
the proof are given in section HI State-dependent QEIs (and related results) have recently 
been studied in an abstract context by one of us [9j, in which they are naturally suggested 
by the mathematical framework. This paper complements that work by giving a concrete 
example of a quantum field which admits a state-dependent bound but cannot admit a 
state-independent one. 

The state-dependent nature of the lower bound raises an important question. It is clear 
that setting 0(f) = — p(f) would provide a rather trivial inequality of the above type. Are 
our bounds similarly trivial? In section [5] we will analyse this question in two ways. The 
first is based on a proposal in [9], in which a QEI would be declared trivial if there exist 
constants c and d such that 

|<p(f)>*|<c + c'|(0(f)>*| (4) 

for all Hadamard We will show that our bound is non-trivial in this sense by considering 
finite temperature states in Minkowski space. The second way uses so-called if-bounds: 
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we show that 12(f) can be bounded by any power of the Hamiltonian greater than 2, while 
p(f) cannot be bounded by powers less than 30 Thus, although the lower bound is state- 
dependent, it is more stringent in energetic terms than any upper bound on the averaged 
energy density: one may say that the 0(f) is of lower order than p(f). In particular, states 
that exhibit large negative energy densities over extended spacetime regions necessarily 
have large positive overall energy. Further comments on the significance of our results are 
given in section [6j 

2 The non- minimally coupled field 

We begin by recalling the definition of the non-minimally coupled scalar field, its quantisa- 
tion and the construction of the stress-energy tensor. This will serve to fix our conventions. 

The classical Lagrangean describing the field on a n-dimensional spacetim^] M = (M, g) 
is given by (CQ), where m, £ are real constants and R is the Ricci scalar with respect to the 
metric g. The constant £ is called the coupling constant. If £ = 0, the field is said to be 
minimally coupled. For £ = (n — 2) /(An — 4) and £ = 1/4, one speaks of conformal and 
super-symmetric coupling, respectively. In this paper we will focus on values £ G [0, 1/4], 
which clearly contains all the special values just mentioned. The Lagrangean (CQ) leads to 
the wave equation 

= 0, (5) 

where := O g + (m 2 + £i?) is the Klein-Gordon operator, and \3 g is the d'Alembertian 
with respect to the metric g. We will follow the standard convention and denote the space of 
compactly supported, smooth, complex- valued functions on M by S)[M). Assuming that 
the spacetime is globally hyperbolic, there is an antisymmetric bi-distribution E^(x,y) 
which is the difference of the advanced and retarded Green functions. 

The theory is quantised by introducing a unital *-algebra 2t^(M), which is generated 
by objects $(/) (/ G &{M)) subject to the relations that (a) the map / — ► $(/) is 
complex linear; (b) $(/)* = $(/); (c) $(P € /) = 0; (d) [$(/), <&{h)\ = iE^f,h)l for 
all /, h G 3>(M). Properties (c) and (d) enforce the field equation and the canonical 
commutation relations respectively. In this framework, states are positive and normalised 
linear functionals (•)#—> C on the algebra 2t^(M). In particular, we will be interested in 
Hadamard states: in such a state the two-point function uf(x,y) = ($(x) ( l ) (y))^ is a 
distribution with a prescribed singularity structure so that the difference between the two- 
point functions of any two Hadamard states is smooth. See [12] and references therein, for 
details on Hadamard states. Some of our later results will be based on a characterisation 
of the Hadamard states in terms of microlocal analysis due to Radzikowski [13]. 

1 The power of 3 emerges by considering a particular family of states, and it may be that p(f) can only 
be bounded by powers of at least 4, as would be natural on dimensional grounds. See section f5. 21 for more 
discussion. 

2 Our sign conventions are those of Birrell and Davies [TD], i.e., the [—,—,—] convention in the classifi- 
cation scheme of Misner, Thorne and Wheeler [TT] . 
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We now turn to the problem of quantising quadratic classical expressions of the form 
G dass (x) = Di((j)®(fi)} c (x), where the Di are linear differential operators on C°°(Mx M) 
with smooth coefficients and [-]c(x) denotes the 'coincidence limit' [F] c (ar) := F(x,x), of 
any smooth function F G C°°(M x M). The quantised normal ordered form of G class in 
the Hadamard state $ is then defined by 

(G*""V*)= [^A:w?:] c (x), (6) 

i 

where :c<jf: = cuf — is the normal ordering of a;* with respect to a reference Hadamard 
state cuq. In Minkowski space one has the distinguished vacuum state Q and therefore one 
usually chooses uj\ = uj^- In all quasi-free representations, normal ordering coincides with 
Wick normal ordering of annihilation and creation operators. 

The classical stress-energy tensor of the non-minimally coupled scalar field can be 
calculated by varying the action of the Lagrangean (TjQ) with respect to the metric, and 
takes the form 

T* ass = (V M 0) (V„0) + \g» v (m 2 2 - (V0) 2 ) + i {g^Dg - V M V, - G, v ) 2 , (7) 

where is the Einstein tensor. The term proportional to the coupling constant £ orig- 
inates in the variation of the coupling term in the Lagrangean density. In a Ricci-flat 
spacetime, the differential operators in ([7]) proportional to £ are still present. So minimal 
coupling and vanishing Ricci scalar result in the same wave equation but in a different 
stress energy tensor. 

To quantise ([7j) in the way we introduced above, we need to bring it into the form used 
in ([6]). This can be done with the definition of the Klein-Gordon operator P%. One finds 
that 

T ciass = (1 _2 O (V M 0)(V^) + ^(l-4O^(mV- W) 2 ) 

- 2£ (V M V,0 + l -R^ 2 - X - (1 - 4£) g, u R4> 2 - g^P^J . (8) 

The last term in the bottom line vanishes "on shell", that is for cf) satisfying the wave 
equation (jSJ). 

We shall often study the energy density of (jSJ) with respect to freely falling observers. 
Assume that 7 is a time-like geodesic parameterised by proper timep i.e., 7 2 = 1 and 
V^y7 = 7 M V M 7 = 0. Using this, together with (jHJ), one can show that the classical energy 
density p^ lass = T^YY on 7 is 

PT S = (1 - 20 (V,0) 2 + ^ (1 - 40 (m 2 2 - (V0) 2 ) 

- 2£ (V 2 + ^R^Ytf ~ \ (1 - 40 R<P 2 - • (9) 



3 We require 7 to be connected, but it does not have to be inextendiblc. 
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Now let T be an open tubular neighbourhood of 7. Take a family of smooth vector 
fields {fj}i=o...n-i on T, whose restriction to 7 is a vielbein with the property that t> | 7 — 7> 
so we have g M ^| 7 = v^Vq — Y^=i v i v i ■ We now introduce the operators 

Pi = ^(V„ ®V„ ) + ^(l-4£) (m\l®l)+^{V Vi ®V Vi )j , (10) 

p 2 = 2(1® S V 2 J, (11) 

p 3 = - (1 ® s R^v v Q )l) + -40(1® SJ R1) + 2(1 ® S P 5 ). (12) 

Here ® s is the symmetrised tensor product, i.e., P ® s P' = {(P(g>P') + (P'®P)}/2. 
Having introduced these operators, one finds that 

pj ass = [p(0®0)] c , (13) 

P = Pi - £P2 + £p3- (14) 

Note that p3(0®0) =0 in Ricci-flat spacetimes if is a solution to the wave equation (JHJ). 
Furthermore, for minimal coupling (£ = 0), we have p = p\. 

We quantise the energy density by replacing the classical point-split field ® by 
the normal ordered two-point function io;* : of some Hadamard state As noted before, 
normal ordering is always performed with respect to some fixed reference Hadamard state 
\EV So the quantised energy density on 7 in the state \I> is simply given by (p quant )^ = 
|j5:a^2 : ] c- Note that our normal ordered energy density is not the same as the renormalised 
energy density (p ren )^ obtained using the Hadamard prescription (see, e.g., [Hj) but they 
are related by (p quan % = (p rcr % - (p ren )* . 

We end this section with a short summary of the non-minimally coupled scalar quantum 
field in the n-dimensional Minkowski space M^j ink = (IR n , rf). We define the measure dp(k) 
by 

d n-1 k 1 

Mk) = (2^pMk)' (15) 

with w(k) = Vk 2 + m 2 and use it to define the one-particle Hilbert space by Ti = 
L 2 (W 1 ' 1 , dp(k)). We will denote the norm and inner product on TC by ||-|| w and (-, -)^, and 
define the bosonic Fock-space J r s ('H) in the usual way. Thus, for each g e 7i, we have an 
annihilation operator a(g) and creation operator aJ(g) with common domain D^N+l) 1 ^ 2 ), 
where N is the number operator on !F s (?i), obeying the canonical commutation relations 
[a(f),a*(g)] = (f, g)nX (recall that g 1— > a(g) is antilinear). 

For any compactly supported distribution / G ^'(M^j ink ), we define /(k) = /(u;(k), k), 
with the Fourier transformation convention 

f(k) = J d n x e ikx f(x). (16) 

If / additionally satisfies the property that \\f\\n < °o and \\f\\n < 00, we can define 

$(/) = a(7) + a f (7) (17) 
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as an operator on D ({N+ I) 1 / 2 ). If we restrict to smooth compactly supported /, the 
operators $(/) restricted to n^L D(N k ) generate a representation of 2l^(M^ ink ) (for any 

o- ' ^ ' 

Formally, we may write 

a(g) = J d/i(k) #(k)a(k) and a\g) = j d//(k) <?(k) a t(k). (18) 

with the a(k) and aJQn) satisfying the commutation relations 

[a(k), a f (k')] = (27T)"- 1 2cu(k)5(k - k')l, (19) 
and the smeared field may be written as 

$(/) = J d n x<S>(x)f(x), (20) 

where 

$(x) = J d/i(k) (a(k)e- lkx + a\k)e lkx ) , (21) 

with x = (t, x) and k = (u;(k),k). Expressions of this type can be made rigorous (cf. 
section X.7 in [13]) and we will make use of this later on. 



3 Quantum states with negative energy density 

In this section we consider the massless quantised scalar quantum field with £ > in 
a (3 + l)-dimensional Minkowski space M^j ink . It will be shown that local averages of the 
energy density are unbounded from below on the class of Hadamard states. 
We start by considering one-particle states for which we have the identity 

:wf:(t,x,lf,x!) = (tf|:$(t,x)$(t / ,x'):¥) = 2Re ((fi|$(t,x)tf)(f2|$(t',x')tf)) , (22) 

as can be shown by writing the fields in terms of annihilation and creation operators. 
The expression for the renormalised energy density becomes relatively simple since the 
differential operator p, defined in (fl4"|) . acting on bi-solutions in Minkowski space is of the 
form 

1 1 3 

P = ® dt>) + ^ (1 " 40 J2( d » ® ^ ~ £ ® 1) + (1 ® %<)) • (23) 

6=1 

So, using (122]) we have 

3 

(*|^™»**)(t,x) = |^(fi|$(t,x)^)| 2 + (l-40^|«9 fe (fi|$(t,x)^)| 2 

6=1 

-4£Re (jSMW^W) x $ (V\ $(t, x) |tf )) . (24) 
6 



For k > 0, we consider (normalised) one-particle states of the form = (v(h K )tt where 
h K (k) = 4ttV2(k - |k| /3)e-^ K /K 2 , for which 

(fi|$(t,x)tt) = J d//(k) e-^l k l- xk ^ K (k). (25) 

Owing to the rapid decay of h K , ty K is Hadamard. It has expected energy 

(V^HVJ = y- (26) 

The functions h K obey h\ K (k) = h K (k/X)/X, as a consequence of which we have the scaling 
relation 

(^ Xk \ p ^^ Xk )(x) = A 4 (^|p quant ^)(Ax). (27) 

Evaluating the energy density of the state at the spatial origin (for example) one finds 
thatS 

<^|p quant * K >(t, 0) = + 8 ^ 2)57r2 {( 3tV + 3tV ) - e(18**« 4 - 44tV + 2)} , (28) 



so, in particular, 

,(2^ 4 

~3tT 

Owing to continuity of the expected energy density and (1211 it follows that to every k > 
there exists a constant r > such 



(vW ua„t^ )(0)0) = _^_^ (2Q) 



(^ip quant ^>(x) < -e^r for a11 x e 5 ( r )' ( 3 °) 

where B{r) is the open ball 

5(r) = {(t,x)|i 2 + |x| 2 <r 2 }. (31) 

Now fix some k and some appropriate r. As a consequence of (j2"7|) we can find a «' > to 
every (arbitrary) r' > 0, such that 

(vD^p^UvXx) < for all x G S(r'). (32) 

In particular, we can take k' = kt/t' . We have therefore constructed a Hadamard state 
with energy density less than — ^K 4 r 4 /(67r 2 (r / ) 4 ) on an arbitrary region B(t') [by trans- 
lational invariance, the same applies to any other spacetime ball of radius t'}. The total 
expected energy of this state is 2kt/(3t'). 

Note that the product of k' and r' is constant. This shows that we may arrange for 
large regions of negative energy density albeit with low magnitude. We may extend the 



''The well-known identity J °° dk k n e k = n\ makes most of the calculations almost trivial. 
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example as follows. Suppose a constant energy density po > is given, and choose an 
integer j > 67r 2 po/(£(/t') 4 )- Then the j-particle state = ty K i ® • • • <8> has energy 
density 

(^\ p ^^)(x) = j(^'|p quant ^'>(^) < -p for all x e B{r'), (33) 
and total energy 

1 K ' ; 3 £(k') 3 1 J £k 3 t 3 ' 1 J 

illustrating that the large negative energy density effects also require large positive overall 
energy, at least in this example. We will see later that this is a general phenomenon. 

Summarising, we have shown that to any bounded subset O of Minkowski space and 
arbitrary constant po > there is a Hadamard state in which the expected energy density 
is less than — p throughout O. In particular any smearing pi uant (f) with a non-negative 
compactly supported distribution f is unbounded from below on the class of Hadamard 
states. 



4 Quantum energy inequalities 

In this section we are going to derive the main result that is to give a lower bound for time- 
like averages of the energy density. We start with a quantum field on a curved spacetime. 
In a second step we specialise these results to Minkowski space, where the vacuum state is 
the preferred reference state. 

4.1 Globally hyperbolic spacetime 

We keep the same assumptions as in section [2j in particular 7 is a time-like, connected 
geodesic parameterised by proper time. Our goal is to find a lower bound for the weighted 
average of the quantum energy density on 7, 

(p quant o 7 >*(f) = J dr f(r) <p quant >*(7(r)) (35) 

in the case where f = f 2 for some real valued function / e f^(R, R). The main task is to 
rewrite (1351) in such a way that the lower bound may be deduced by discarding manifestly 
positive terms. 

To start, we use (fill) to get 

(P quant °7>*(f) 

= ([p! *,»:], o 7) (f) - £ ([p 2 :u,*:] c o 7) (f) + £ ([p 3 :^:) c o 7) (/ 2 ). (36) 

Each term on the right-hand side will be treated in turn. It will also be useful to de- 
fine (p(r,r r ) = {j(t), 7(t')), and to write <p*F to denote the pull-back ip*F(r,r') = 
F(j(t),j(t')) of a smooth function F from M x M to R x R. 
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The first term on the right-hand side in (|36|) may then be rewritten, following [5], as 






drdr' e^^f{T)f{r')^{ Pl :^:)(r,r') 



(37) 



where f a (r) = e lOLT /(t). Here, we have made use of the Fourier representation of the 
5-function and also the symmetry of the normal ordered two-point function to arrange 
that the a-integral takes place over the positive half-axis. A decomposition of this kind 
will be referred to as a point- splitting trick, see [5]. The distributional pull-backs of the 
form v ? *(pi ^f) appearing in the last step were shown to exist in [5j, using the microlocal 
characterisation of Hadamard states given in [13"] . Moreover, if £ < 1/4 then these dis- 
tributions are positive type, i.e., <p*(pi ujf)(C,C) > f° r an C ^QR). This is a direct 
consequence of Theorem 2.2 in [5] and the form of p\. In particular, each integrand in 
the bottom line in (13~T|) is non-negative. Finally, the integrals converge, because (as shown 
in [5]) (p*(pu2){f a , fa) is of rapid decay asa-+ +oo for any Hadamard state ^ and partial 
differential operator p with smooth coefficients. 

To treat the second term on the right-hand side in ( 13"6l) . we will need the following 
identity, proved in appendix [At 

Theorem 4.1. Let F be a smooth function on M x M and d be a partial differential 
operator of the form d = C M V M; with a smooth vector field (. Then 



for h G @(M,R). 

Now choose a function f? G M) such that /r ° 7 = /■ Applying theorem 14.11 

with F = luf '•, h = fq- and d = C^V^, where ( is some smooth vector field on M with the 
property that £| 7 = 7, yields the identity 



2/2 [(1 0S V?):w*:] c + V c [(l ® s V c /|):^: - (1 ® s #V f ):<«#:] c 
= -2[(V f / r ® s V f / r ):wMc + 2V c /|[:a;*:] c + 2V C [(/ T ® s V c /r):^|:] c . (39) 



2/i 2 [(l ® s <9 2 )F] C + 9(1®, <9/i 2 )F - (1 ® s /i 2 <9)F 



J c 



= -2[{dh ® s a/i)F] c + 2(<9/i) 2 [(l ® s 1)F] C + 29[(/i ® s dh)F] 



CI 



(38) 
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All expressions are well defined, since :u;*: is smooth. The first expression on the left-hand 
side is nothing but /f-f/^:^* :] c . Let us turn to the other terms. Recalling that r is the 
proper time along 7, we can write 

J dr (V4(l ® s V^):u;*:] c o 7) (r) = J dr d T ([(1 ® s V^|):^|:] c o 7) (r), (40) 

which vanishes, since fr is of compact support. Further terms in (|39|) will vanish for the 
same reasons after integration. We finally obtain that 



dr/ 2 (r) ([p 2 :^:] c o 7 )(r) = 2 J dr (d T f) 2 ^:w»:(r, r) 

-2 /" dr ([(V^/T® s V^/r):a;*:] c o 7 )(T). (41) 

The second integral on the right-hand side may be rewritten (up to a factor), using the 
point-splitting trick, as 

J dr ([(V^/ T ® s V^T):^:] c o 7 )(r) 
= J" 7« 2 / d ^r' Ta{r)Ur') ^:w*:(t, r') 

= r /a ) _ r .^qt /a)) (42) 

io 1 Jo 71 

where we have also used integration by parts in r, r' and the fact that / is of compact 
support. As before, the bottom line of (142]) is a difference of two non-negative terms. 

Finally let us put the results of (1371) . (j4"Tj) and (l4"2l together with the remaining term 
in (IMjl FI We find that 

(p quant °7>*(/ 2 ) 

r 00 d a f°° da 

= / — <P*(Pi u 2 ) (ja, fa) - / — <P*(PlU%)(f a ,f a ) 
Jo 7T io ^ 

/* OO 1 /* OO 1 

+2 £ / ^a 2 / ) - 2£ / —a 2 ^wg(^, / a ) 

-2£ I dr (d T f) 2 ¥»*:a;»:(r, r) + £ | dr f(r) <p* (r, r), (43) 

which, noting that <p*:a;f:(T, r) = (:$ 2 : o 7)^(r), results in the following 

5 We do not apply the point-splitting trick to this term as we cannot necessarily find smooth real square 
roots of the geometrical quantities involved. 
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Theorem 4.2. Let uj® be the two-point function of a reference Hadamard state for the 
non-minimally coupled scalar field with coupling constant £ £ [0,1/4], defined on a glob- 
ally hyperbolic spacetime with smooth metric. Furthermore, let 7 be a time-like geodesic 
parametrised in proper time r and let f G i^(R, R). On the set of Hadamard states, we 
then find 

(p«-o 7 )(/ 2 )>-^(/), (44) 

where 

&{f) = Qi(/)1 + £ (:$ 2 : o 7) (0 B [/]) + £ (:$ 2 : o 7) (££[/]), (45) 

= r — ^ 2 °) (7a, fa) + 1t0? V*U° 2 (U, f a )] , (46) 

JO 7r 

and 0.3(f) and £r c (f) are functions in f^(R, R) given by 

Q B [f](r) = 2{d T f{r))\ (47) 

&M(t) = f(r) 2 (^7V-^(l-40i?)(r). (48) 

Furthermore Q. A (f) and £2b[/] are non-negative. 

This result follows from the previous discussion by discarding manifestly positive terms. 
We remark that £l^ A (f) depends on the reference state and that for £ = 0, we recover results 
known for minimal coupling [5]. Moreover, H^(/) vanishes if the region of interest is Ricci- 
flat. 



4.2 Minkowski space 

In this section we apply the results derived in the previous subsection to n- dimensional 
Minkowski space. Without loss of generality, we average in the time argument r = t at the 
spatial origin, i.e., j(t) = (r, x ). We choose our reference state to be the vacuum state 
Q, which has two-point function 

^(t,x,t',x') = J d/i(k)e- i[( '-'' ) " (k) - (x - x ' )kl , (49) 

in the distributional sense. For g £ ^(R), we find that 

1 <? r°° k n ~ 2 

where S , „_ 2 is the surface area of the (n — 2) dimensional standard unit sphered We also 
have the identity 

n-1 

m 2 (p*u$(g ®g) + J2v* (($ ® diV" ) {g ® g) = ^ ((<9 ® ® <?), (51) 

i=l 

6 We have S^n = 2v / 7r™/r(m/2), where T is the Gamma function. 
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which follows from the spacetime translation invariance of the vacuum and the field equa- 
tion ((Sj). So one can absorb the mass term and the spatial derivatives appearing in the 
definition of Q. A (f) (via pi) into a further term that involves time derivatives. One finds 
that (HU) becomes 



Theorem 4.3. For the non-minimally coupled scalar quantum field in n- dimensional 
Minkowski space M^ injfc , 

(p^o 7 )(/ 2 ) > = - (5i(/)l + £(:$ 2 :° 7 ) (Ob[/])) , (52) 

in the sense of quadratic forms on Hadamard states, where 

& A (f) = J^J da J dk ^) (C 1 " 2 £V 2 W + 2^ 2 ) /(a + (53) 

and 

n B [/](t) = 2(^/(t)) 2 (54) 

/or / G 0(R, R) and £ G [0, 1/4]. 

It is easy to see that 0. A (f) is non-negative for £ G [0, 1/4] and that 0^ _o (/) = Q.(f), 
where 0.(f) is the lower bound found in [16] for the minimally coupled (£ = 0) scalar field. 
As a consequence, theorem 14.31 recovers the results of [16] for minimal coupling. 

We can write jQ^Cf) in the fornfl 



q poo 



X ( -Q n ,2 (—) ~ ^-^—Qn,l (—) + ^^—Q n , (—)), ( 55 ) 



n " ' \m/ J n — 1 ' \m/ n — 2 \;// 
where the non-negative functions Q n .fc are defined by 

QnM = [ dx (x 2 - 1)M/V, (56) 

for n + > 2 . They vanish for k + n = 2 and else have the properties that Q n .fc(l) = and 
Qn,k(y) 1 as V - °° • As they are continuous, they are therefore bounded. It might be 
useful to use the following estimate for £l A (f), which follows from the previous discussion 
for n > 2 and £ G [0,1/4], 

&!(/) < /"^ |/| 3 (uK- (57) 

AW ; ~ (27r) m 2n(n - 2) J v ; 

This estimate is true for the massive and the massless case and one can find a similar 
estimate for the two-dimensional case. 



7 We will assume that n > 2 and m > 0, but one can find similar expressions for these cases as well. 
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To conclude this subsection, we investigate the behaviour under rescaling of the aver- 
aging function. The smearing function f\(t) = f(t/X)/y/\, for A > 0, has the property 
that || fx || L i = || /|| £2, and its Fourier transform satisfies the identity 

(A(«)) 2 = A(/(A«)) 2 . (58) 

One can conclude from this that 0^(/a) = 0(A~ n ) as A — > oo. In fact, even faster decay 
could be concluded if m > using the arguments of [I7j . It follows that \£r A (f\) — > as 
A — ► oo. For states with |(:$ 2 : o j)^(t)\ < c(l + |t|) 1_£ , for some positive constants c, e one 
can show that A(:$ 2 : o 7)#(0.b(/a)) — ► as A — > oo, and we therefore obtain the averaged 
weak energy condition (AWEC) for £ 6 [0, 1/4] in the form 

liminf ( dt (p quant o 7 ) v (t) /(t/A) 2 > 0. (59) 

This is in line with a result in [18], which shows that AWEC holds for states in which the 
particle number and the energy is bounded (see the penultimate paragraph in section III 
of [E]). In the minimally coupled case, AWEC is also known to follow from QEIs [T9l fT6]. 
We will return to the AWEC briefly in section 15.21 



5 Investigation of state-dependence 

The lower bounds that we derived have the characteristic that they are state dependent 
except for minimal coupling. As all previously known QEIs are state independent it is 
important to understand the nature of the state dependence to ensure that our bounds are 
not vacuous. 



5.1 KMS states and temperature scaling 

In this subsection we analyse the temperature scaling behaviour of the stress energy tensor 
and the bound in theorem 14.31 for a KMS state \f^, i.e., a thermal equilibrium state at 
positive temperature as seen by the observer on 7. Its two-point function uj^ in a 
n-dimensional Minkowski space, with n > 3, is given by 

f /e -i((t-t>(k)-(x-x')k) e +i((t-t')w(k)-(x-x')k) 

J 2 (t, x, t>, x') = y d/x(k) ( i - e _^ (k) + ) . (60) 

We renormalise the two-point function of the KMS-state by subtracting the two-point 
function of the vacuum (T49|) . In the coincidence limit, we find that 



:u%:] c (t,x) = / d/i(k) 



e /Mk) _ 1 

B n>0 {/3m), (61) 
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with the positive function B n ^ defined on [0, oo) for k > by 

Q POO 

Bn Ja) = ^- I dz (z 2 -a 2 )^ 



(27r) n - 

The expression ( IBTj) is positive and invariant under spacetime translations. Thus the state- 
dependent part of the lower bound in theorem 14.31 is given by 

(:$ 2 : o tV(jQ b [/]) = 2p 2 ' n B nfi (Pm) \\ff L2 , (62) 

while the state independent part 0^(/), is obviously independent of (3. On the other hand, 
the renormalised energy density of this state is 



d/*(k)- 



(p^ uant Mt,x) = [p:^:] c (t,x) 

2a; 2 (k) 

= (3~ n B n ^m), (63) 
so that the time averaged energy density is given by 

(p quant o 7 M/ 2 ) = /r n £n, 2 (/3m) ||/|| 2 L2 • (64) 
We can now state the non-triviality result: 

Theorem 5.1. The bound for the energy density of a non-minimally coupled scalar quan- 
tum field given in theorem \4-3\ is non-trivial in the sense of ' t 9], i.e, there do not exist 
constants c, d such that 

|(p 9 " ant °7>*(/ 2 )| <c + c'\{&(f)}*\ (65) 
for all Hadamard states ^ unless f is identically zero. 

Proof. From the previous discussion, we find that for a fixed non-trivial smearing func- 
tion /, in the limit of high temperatures 

lim/5"(p^ ant o 7 )^(/ 2 ) = B ni2 (0) \\f\\ 2 L2 > 0, (66) 
and lim/rniCf) = /T(:$ 2 : o 7 }^(Ob[/]) = 0. (67) 

Now if the bound in theorem 14.31 is trivial, then there exists some constant c such that (j4]) 
holds. This implies that 

< B n>2 (0) lUfe < hm/T (c + d \& A (f)\+d |£(:$ 2 : ° 7>*/»(O fl [/])|) = 0, (68) 

which is a contradiction. □ 

A more refined formulation of this result would be that our bound is non-trivial with 
respect to high temperature scaling for KMS states. 
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5.2 Energy behaviour 

The previous results already show that the lower bound in theorem 14.31 has a different 
scaling behaviour from the energy density itself. We will now present a more general 
analysis that gives more insight into this, again working in four-dimensional Minkowski 
space M^ ink . 

The lower bound in theorem 14.31 only depends on the fields :$ 2 : and 1. On the other 
hand, the energy density also involves terms such as :$ 2 :. This results in a crucial difference 
in their energy behaviour. To be more precise, we seek values of p, q e M. + for which there 
are constants Cf >q ,df such that 

c,, a (H + ml) q > (p quant o 7 )(/ 2 ) > -c' f}P (H + mlf (69) 

holds (in the sense of quadratic forms) on the set of Hadamard vector states. In the 
minimally coupled case we know that there is a state-independent lower bound so we 
may take p = 0; however for £ e (0, 1/4] we have already shown that (p c i uant o 7)(/ 2 ) is 
unbounded from below, so p must be strictly positive if ( 1691) is to hold. By theorem 14.31 it 
is enough to show that < c 'f p (H + mt) p to conclude that the right-hand inequality 

in (169]) holds; we will show that this is possible for any p > 2. 

On the other hand, we will show that the left-hand inequality in (|69|) cannot be satisfied 
for q < 3. In this sense, our lower bound represents a non-trivial constraint. Indeed, 
the situation here is reminiscent of the sharp Garding inequalities studied in the theory 
of pseudodifferential operators, in which operators with positive classical symbols may 
be bounded from below 'with a gain in derivatives', i.e., by an operator of lower order. 
Although the analogy is not direct, it seems worthy of further investigation. 

As a consequence of this analysis we immediately obtain another proof of theorem 15.11 
namely that the lower bound in theorem 14.31 is non-trivial in the sense of [9] (but this time 
using states in the domain of a power of the Hamiltonian rather than KMS states). 

We begin by establishing our claim relating to the left-hand inequality in ( 1691) . To do 
this we first consider the massless field and construct the following one-particle state 

m K = ^L J d/i(k)e~ |k|/K a t (k) a (70) 
It is straightforward to calculate that (iiP)# re = (j + 1)! and tha 



^■«=?|iw-Wr (71) 

We have a peak at t = 0, where the expectation value scales like the fourth power in the 
Hamiltonian. However, this pointwise behaviour will not hold for the smeared field. Let 
us assume that / is an non-trivial, integrable function. We find that 

Um = lim ^(2 + 3{)A = 2(2 + 30 (72) 



8 For the concrete calculation in (fTTj) we assumed for simplicity that the geodesic 7 is located at the 
spatial origin. 
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i.e., we have an asymptotic scaling like the expectation value of H 3 . In the limit where k 
becomes large, the high momenta are dominant. If the field is massive, the mass would 
therefore become negligible. We can deduce that the same scaling behaviour remains 
true for the massive field, at least asymptotically. So the smeared energy density for the 
scalar field in Minkowski space scales asymptotically at least with the third power of the 
Hamiltonian, i.e., the left-hand inequality in fl69l) can only hold if q > 3. It is worth 
remarking that we could have conducted the same analysis using the one-particle states 
investigated in section [3], and with the same result: although the energy density in these 
states is negative at the spacetime origin, one may check that the smeared energy density 
is positive for sufficiently large k. 

The second part of the discussion aims to establish suitable i7-bounds on the state- 
dependent part of our lower bound, i.e., the field :$ 2 :. if-bounds have been discussed 
elsewhere (see e.g., |20j), but we require more detailed information on the power of the 
Hamiltonian involved and on the controlling constants than we have been able to locate in 
the literature. The following discussion may therefore be of independent interest. 

We will now assume that m > 0, but allow general spacetime dimensions n > 2. Let h 
be the one-particle Hamiltonian and dr be the second quantisation map (so, for example 
the Hamiltonian is H = dT(h)). Let us initially restrict to the domain Dy C T S ^H), 
defined as the space of vectors in Fock space all of whose n-particle wavefunctions are 
Schwartz functions and all but finitely many of which vanish identically, see section X.7 
in [15] . This is a dense domain in the Fock space and for every vector \& G Dy, we find 
that k ^ ||a(k)#|| 2 e ^(M™" 1 ) and that 

||dr(/i p ) 1/2 ^|| 2 = J d/i(kK(k) l|a(k)*|| 2 (73) 

for any p G R. 

Let g £ yiW 1 ^ 1 ) C H. Then uj~ p / 2 g E H, where u acts on H by multiplication and 
the Cauchy-Schwarz inequality implies that 

\\a(g)y\\ 2 < y'd^(k)d^(k , )|<a(k)*|o(k , )*>lb(k)|b(kOI 

< d/i(k)||a(k)*|||<7(k)| 

J d/i(k) (^ 2 (k) \\a(k)n)\^- p/2 g(k) 

< ||dr(/ i p ) 1/2 ^|| 2 - \\uj- p/2 g\\ 2 H pel, (74) 

for ^ E Dy. Due to the construction of Dy, one can find that Dy C D(H P ) for any 
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p 



E R. Now let v&W C Dy be a /-particle state. We have 



(H p ^) (k u . . . ,k L ) = (X>( k *)J * (0 (ki,...,k z ), 

i 

(dr(^)^)(k 1; ...,kO = ^^(ki^Ckx,...,^). 

i=l 

For a, 6 > and p > 1 we find that (a + 6) p > a p + &p@ so 

' I \p i 

$>(k.) >$>^). 



(75) 
(76) 



.1=1 



i=l 



This implies that < dr(/i p ) < H p on Dy for p > 1. Using this, the commutation relations 
and (|74l) . we recover if -bounds of the form 



\\a(g)Vf < \\H p/2 ^\ 
||a%)^|| 2 < \\H p / 2 ^\ 

for p > 1 on Dy. 

Now define the distribution space, 



\"- p/2 9\\m 

\"- p/2 9\\n + \\9\\n\\n 2 , 



(79) 



C q = {F G ^(M- ink )|||F||^ := || W *F|& < oo}, (80) 

where ||-|| is a semi-norm. Since the field is massive we have the inclusion £ C £ g for 
g < 0. The field ^(i 71 ) defines an operator on Fock space on the domain D ((N + l) 1 ^ 2 ), if 

F,Fe£ (or equivalents F,F G H). Since Dy C D ({N + l) 1 / 2 ), $(F) is well-defined 
on Dy. 

Now assume p > 1, ^ G -Dy and F,Fe Cq. Using the inequality (w + 1>) 2 < 2(u 2 + t> 2 ) 
we find 

||$(F)^|| 2 - ||$(F)n|| 2 



< 



a(F)V + a f (F)^ 
< 2||F p/2 ^|| 2 - fllFll 2 ln + \\F\ 

— II II \\\ M-p/2 11 1 



2 

-p/2 



+ ll^lloll^ll 2 



Since this inequality is valid on Dy, which was dense in the Fock space !F S (TC), this result 
extends to all * G D(H p / 2 ). 



9 Since the function x — > x p 1 is monotone increasing for p > 1, we have 

(a + b) p = a(a + b)^ 1 + b(a + bf' 1 > a? + W. 



(77) 
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In order to apply these results in conjunction with the point-splitting trick, we have to 
establish a connection between pulled back two-point functions and their representation 
in terms of distributionally smeared fields acting as operators on a certain domain. In 
particular one has to check that for some F(x) = ((t) ® <5 X0 (x) with ( G f^(R), we have 

i^)HO = \MF)n 2 . (82) 

The left-hand side is well-defined as discussed before. The right-hand side is well-defined 
as F, F G Co. The identity can then be shown by constructing a sequence F r G ^(M^ ink ), 

with u%(F r ,F r ) -> ((p*uf )(£,() and F\. —> F = C®S> in ^0 The latter property 
ensures that $(F r )\l/ — > $(£ ® 5 XQ )\I/ and the identity therefore holds because uf(F r , F r ) = 
||$(F r )^|| 2 for test functions F r G ^(M^ ink ). 

We are now able to apply the above result to find if-bounds on the Wick square, 
smeared along the inertial curve 7. Applying the point-splitting trick and ([8*1]) . we find 
that 

<:$ 2 :°7>*(/ 2 ) 

= — (\Mfa ® <U*f " \Wa ® SjOf) 

JO ^ 

< 2|i^/ 2 vi/|| 2 . r~ (ii^®^ii 2 _ p/2 +ii/ a ®5 X0 ii 2 _ p/2 ) 

»/ 



, r°° da 9 
+ 11*11" / — ll/<*®M 
Jo * 

2\\ H ^\\ 2 . r-\\f a ®s X{ f_ p/2 +un 2 r-n/«®us 

J-00 7i JO 71 



(83) 



where / G f^(R, R) and we used fl82|) with F(x) = fit)e tat ® 5 Xo (x). For convenience, let 
us introduce the positive quadratic functionals 

5- P / 2 (/) = 2 ^ ^ ||/a ® <U 2 _ p/2 (84) 

Co(/) = f — ll/a®«Ui!- (85) 
Jo 7r 

Lemma 5.2. Let (p + 2) > n ; where n is the spacetime dimension, and let f G J^(R, R). 
For £/ie massive case, then Coif) < 00 anc ^ ^-p/2(/) < 00 ■ 

Proof. To show that |Co(/)| < 00, we see that 

/r®X (k) 2 = />(k) + a) 2 < ° +1 < * w _° , (86) 



'One can do this by defining F r (t,x) — £(i)xr( x )j with the approximate identity \r € @(M. n x ). 
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for some positive constant c, where we made use of the fact that / is in Schwartz space. 
It follows that there exists a constant c' such that 



Ilia® <5xo Ho < 



c 



m + a 



\2-' 



17) 



which is integrable in a on IR + proving that Co(/) < oo. 

To show that LB_ p / 2 (/) < oo, realise that there exists a constant c, such that 



< 



[w? + (w(k) + a) 2 )' 



(88) 



Using this inequality, we get 
da 



7T 



ll/a®5x ||_p/2 



< 



d/i(k) CU' P (k) 



da 



7r (m 2 + (a?(k) + a) 5 



< - / d/i(k) ^(k), 
m J 



19) 



where we made use of Tonelli's theorem. The last expression, however, is finite due to the 
restriction on p. □ 

We note that, under the assumptions of lemma l5\2l the quantities in ( 184"1) and ( |85j) may 
be estimated by similar arguments to those used to obtain (157|) . Combining inequality (185]) 
and lemma I5.2[ we can now summarise our result as the following if -bound: 

Theorem 5.3. Let n be the spacetime dimension, p > (n — 2), / G ^(R, R) and ^ G 
D{Hp' 2 ). Then 



(:$ 2 : o 7 >*(/ 2 ) < B-j/aCf) ||iP /2 vI>|| 2 + C (/) ||*f < oo. 
As an immediate consequence of theorem 15.31 we have: 



(90) 



Corollary 5.4. Subject to the assumptions and notation of theorem \4-3\ and theorem \5.3[ 

we have 

\(&%(f)\ < + 2£C ($/)) ||vl/|| 2 + 2^ p / 2 (9 t /)||^/ 2 vl/|| 2 . (91) 

It follows that there exists a constant c'j p for which 

(p quant o 7 )(/ 2 ) > 5 (/) > -c' f>p (H + miy (92) 
as an inequality of quadratic forms for Hadamard vector states 
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In particular, we see that large negative time-averaged energy densities (for a given 
smearing function) can only be obtained at large positive overall energies. 

As an application of these results we may again consider the rescaled test function f\ 
as in our earlier discussion of the AWEC. One can show that \Co(dtf\) and \B- p /2(dtf\) 
both tend to zero as A —>■ oo (for the allowed values of p) which, together with our earlier 
observation that AH^(/a) — > 0, yields a proof of AWEC in the form (159|) for all Hadamard 
vector states (which necessarily belong to the domain of H p / 2 ). 

Now let us return to the four- dimensional case. Our result shows that a bound of the 
type fl92l) can be satisfied for any p > 2 (although we cannot exclude the possibility that 
it might also hold for some p < 2). The important point is that our QEI is a non-trivial 
restriction on the averaged energy density because the upper bound in fl69l) cannot be 
satisfied for any q < 3 (although it can be satisfied for any q > 4 by adapting our if-bound 
arguments) . 

6 Conclusion 

Our main results may be summarised as follows. First, we have shown (at least for massless 
fields in four-dimensional Minkowski space) that the non-minimally coupled scalar field 
with £ > admits states with arbitrarily large negative energy density over arbitrarily 
large bounded spacetime regions. Thus the non-minimally coupled field cannot obey QEIs 
of the type previously studied in the literature, in which there is a state-independent 
lower bound. Second, by combining the QEI derivation for the minimally coupled field [5J 
with techniques previously applied to the classical non-minimally coupled field [3], we 
have derived a new type of QEI with a state-dependent lower bound, for couplings in the 
range (0, 1/4] in general globally hyperbolic spacetimes. Third, we have analysed the state- 
dependence of the bound, which involves the Wick square of the field, rather than the Wick 
powers of its derivatives that appear in the energy density. This involved the formulation 
of various ^/-bounds, which may be of independent interest. Roughly speaking, the results 
of this analysis tell us that the lower bound scales more softly with the overall energy scale 
than the energy density itself. Thus we may conclude that negative energy effects with 
large magnitude, while possible over large regions, require more energy to achieve than 
positive energy densities of the same magnitude and that the 'energy budget' for these two 
effects will grow with a different power. 

In the light of our results, it seems reasonable to expect that generic interacting quan- 
tum fields will not obey state-independent QEIs, as has also been argued on physical 
grounds for a particular model in [21]. The state- independent QEIs that have been found 
for other free fields and conformal fields in two dimensions (see [22J) should therefore be 
regarded as particularly simple cases. In general, then, the aim should be to establish 
non-trivial state-dependent bounds. Suitable generalisations of our if-bounds to higher 
order Wick polynomials may be useful in this context. 

It also becomes important to understand whether one can draw useful physical con- 
clusions from state-dependent QEIs. For example, the state-independent QEIs have been 
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used to place constraints on exotic spacetimes [231 1211 [25j [26] and are linked to questions 
of thermodynamic stability [271 12H] • A first indication that state-dependent bounds can be 
used for similar purposes may be found in our derivation of AWEC. In general we expect 
that similar results will obtain for state-dependent bounds, once one restricts to states 
of a given energy scale. Likewise, we would expect the phenomenon of quantum inter- 
est [29l [30] to be governed by the energy scale in this context. Finally, our general QEI 
is of the so-called 'difference' type: it constrains the normal ordered energy density rather 
than the (Hadamard) renormalised version. This restriction has recently been removed for 
the minimally coupled field in [31] to obtain an 'absolute' QEI; one would expect that this 
can also be adapted to the non-minimally coupled field. 



A Proof of theorem 4.1 



Throughout the appendix [-] c denotes the 'coincidence limit', i.e., the diagonal of smooth 
functions on M x M, but for the sake of clearity we will omit the arguments. 

Proof of theorem 4-1 First, we look at the top line in ( |38l) . To make the calculations 
clearer, we start by calculating the expression without the symmetric product. We know 
by Synge's theorem that 

d[F] c =[(d(g)l)F] c +[(l®d)F] c , (93) 

where F is a smooth function on M x M, so 

d[(t <g> dh 2 )F - (1 <g> h 2 d)F} c 
= [(d <g> dh 2 )F + (1 ® d 2 h 2 )F] c - [{d ® h 2 d)F + (1 ® dh 2 d)F] c 
= l(d®(dh 2 ))F+(d®h 2 d)F] c 

+ [(1 ® (d 2 h 2 ))F + 2(1 ® (dh 2 )d)F + (1 <g> h 2 d 2 )F] c 

-[(<9 (8) h 2 d)F} c - [(1 ® {dh 2 )d)F + (1 <g> h 2 d 2 )F} c 
= 2(dh 2 )[(l® s d)F} c + (d 2 h 2 )l(l® t)F] c . (94) 

Where we want to emphasise that there is a symmetric product in the last line. Now it 
is quite obvious that the same calculation is valid if one exchanges the arguments of F. 
Thus adding the term 2h 2 [{l ® s d 2 )F] c on both sides, we find the following identity for the 
symmetrised expression 

2h 2 [{l ® s d 2 )F} c + d[(l ® s dh 2 )F - (1 ® s h 2 d)F} c 
= 2h 2 [(l ® s d 2 )F} c + 2(dh 2 )[(l ® s d)F] c + (d 2 h 2 ) [(1 ® s 1)F] C . (95) 

Here the top line is obviously identical with the expression in the top line of (|38|) . 
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Now let us do the analogous calculations for the bottom line in (!38|) . As before, we 
start for clarity without the symmetric product 

2d[(h <g> dh)F] c - 2[(dh <g> dh)F] c + 2(dh) 2 [(l <g 1)F] C 
= 2[(dh <g> + (/i ® <9 2 /i)F] c - 2[(d/i ® 5/i)F] c + 2(«9/i) 2 [(l ® 1)F] C 
= 2h[(l ® (<9 2 /i))F + 2(1 ® {dh)ff)F + (1 ® /i<9 2 )F] c + 2(<9/i) 2 [(l <g> 1)F] C 
= 2/i 2 [(l <g> <9 2 )F] C + 2(2/i<9/i) [(1 ® 9)F] C + (2/i<9 2 /i + 2(9/i) 2 )[(l ® 1)F] C 
= 2h 2 [(l®d 2 )F} c + 2(dh 2 )l(l®d)F} c + (d 2 h 2 )l(l® 1)F] C . (96) 

Again, if we symmetrise this, we get 

2d[(h O s dh)F] c - 2[(dh ® s dh)F] c + 2(dh) 2 [(l ® s 1)F] C 
= 2/ i 2 [(l® s 9 2 )F] c + 2(9/i 2 )[(l® s 9)F] c + (9 2 /i 2 ) [(1® S 1)F] C . (97) 

Comparing the bottom row of fl97|) with the bottom row of fl95l) proves the theorem. □ 
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